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I.  INTRODUCTION 

A 

This  is  the  final  report  on  a  1-year  study  of  adaptive  arrays  in 
a  scattered  jamming  environment.  It  is  assumed  that  the  jamming  signals 
are  scattered  into  the  main  receiving  beam  of  a  radar  or  communication 
system.  The  scattering  medium  could  be  terrain,  rain,  or  chaff.  Con¬ 
ventional  adaptive  nulling  which  uses  angular  discrimination  cannot 

be  used  against  this  type  of  interference  without  also  nulling  the  desired 
signal,  since  both  the  scatterers  and  signal  source  are  in  the  main  beam. 

A  method  of  cancelling  the  scattered  jamming  signals  is  outlined  in  See-;~2 


this  report.  The  theory  of  this  type  of  adaptive  cancellation,  presented 
in  earlier  reports  on  this  contract,  is  also  reviewed, -in  Sec»-2— -  r 

The  technique  for  cancelling  scattered  jamming  uses  delayed  replicas 
of  the  jamming  signal.  Since  the  area  or  volume  of  scatterers  may  be  large, 
the  resulting  adaptive  canceller  may  require  a  large  number  of  adaptive 
weights  (e.g.,  1000  or  more).  This  suggests  a  complicated  adaptive  system 
and  the  importance  of  devising  methods  of  simplifying  the  adaptive  weight 
computation.  It  has  been  shown  that  the  required  number  of  adaptive  degrees 
of  freedom  depends  on  the  frequency  spectrum  of  the  jamming.  With  wide 
band  noise  jamming,  this  frequency  spectrum  depends  on  the  receivers  which 
can  be  designed  to  facilitate  scatter  cancellation.  In  earlier  reports, 
two  frequency  spectra  have  been  investigated,  the  rectangular  bandlimited 


spectrum  and  the  Gaussian  spectrum.^These  results  are  reviewed  in  Secs.  3 


and  4  of  this  report.  During  the  last  quarter,  a  bandlimited  cosine 
spectrum  was  investigated.  These  recent  results  are  contained  in  Sec.  5. 

It  is  shown  that  the  cosine  spectrum  requires  fewer  adaptive  degrees  of 
freedom,  so  that  simpler  implementations  can  provide  20  to  30  dB  of  scatter 
cancellation. 
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During  this  study,  the  case  of  two  jammers  was  also  analyzed.  It 
was  shown  that  the  same  technique  for  scatter  cancellation  can  be  used 
against  two  jammers  when  both  are  illuminating  a  scattering  medium  in  the 
main  beam.  Results  for  the  two  jammer  cases  are  contained  in  the  first 
quarterly  report  on  this  study 


2.  BACKGROUND 


Adaptive  array  antennas  can  provide  a  significant  improvement 
in  the  performance  of  communication  and  radar  systems  when  jamming  is 
present.  The  weights  on  the  array  elements  are  adaptively  controlled 
to  place  antenna  pattern  nulls  in  the  directions  of  jammers,  thus 
improving  the  system  signal -to-noise  ratio.  This  adaptive  array 
technology  has  been  experimentally  verified  and  is  currently  being 
incorporated  in  many  radar  and  communication  systems. 

When  a  system  utilizes  a  low  sidelobe  receiving  antenna  plus 
adaptive  nulling,  a  high  degree  of  immunity  to  direct  line-of-sight 
jamming  is  achieved.  In  these  cases,  a  second  type  of  interference 
will  sometimes  limit  the  system  performance,  viz,  jamming  scattered 
from  terrain  or  chaff  into  the  main  receiving  beam.  This  interference 
is  received  from  the  desired  look  direction  and  cannot  be  rejected  by 
angular  nulling. 

A  method  of  cancelling  scattered  jamming,  which  has  been 
described  in  detail  in  earlier  ASI  reports,  is  summarized  in  Fig.  1. 

A  jammer  is  illuminating  a  scattering  area  or  volume  in  the  main 
beam  of  a  receiving  antenna.  The  scattered  jamming  is  delayed  by 
(R.j+R2)/c,  which  extends  over  a  time  interval  depending  on  the  geometry 
as  illustrated  in  Fig.  1.  This  same  jamming  signal,  delayed  by  Rd/c, 
is  received  by  an  auxiliary  antenna  element  located  near  the  receiving 
antenna.  A  delay  line  is  used  to  equalize  the  two  delay  paths.  When 
the  scattering  region  is  extended  in  area,  an  interval  of  delay  is 
covered  by  taps  on  the  delay  line.  These  taps  correspond  to  an 
interval  bracketing  the  interval  of  (R-|+R2-Rd)/c ,  and  are  spaced  by 


Fiaure  1  —  Scatter  Cancellation 


roughly  1/B,  where  B  is  the  jammer  bandwidth.  The  weights  w-j  ,w2 , . . .  ,wN 
are  selected  to  minimize  the  jamming  residue  in  the  output. 

As  in  Fig.  1,  let  y  denote  the  output  of  the  main  antenna  and 
z  the  final  output  after  cancellation  of  the  scattered  interference. 

The  tapped  delay  line  outputs  are  denoted  by  the  column  vector, 

E  =  {Ep  E2,  ....  EN)T.  Let  W  =  •••»  wN)T,  denote  the  set 

of  adaptive  weights  applied  to  the  delay  line  taps.  Then, 

Z  =  y  -  W*E  (1] 

It  has  been  shown  that  the  weights  which  minimize  the  scattered  jammer 
power  in  the  output  z  are 


W  =  M"1  Ey*  (2! 

* 

where:  M  =  EE  =  covariance  matrix  of  tap  outputs 

Ey  =  column  vector  with  component  Eny* 

With  the  optimum  weights  of  (2),  the  output  power  is 

|Z1I=  \y\2  -  E*y  M'1  Ey*  (3] 

Let  p(t)  denote  the  autocorrelation  function  of  the  interference 

signal . 


p(t)  =  f(t)f*(t-t) 


(4) 
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Then,  the  covariance  matrix  of  tap  outputs  is 


p(0)  p(A)  p(2a) 


p((N-1)A) 


P ( A )  p(  )  p(A) 


p((N-1)a) 


where  A  is  the  spacing  between  taps.  This  matrix  is  Toeplitz  and 
Hermitian.  In  most  cases,  the  spectrum  is  symmetric  about  the  carrier 
frequency  so  M  is  real  and  symmetric.  The  Toeplitz  property  [i.e., 
that  Mmn  is  a  function  only  of  (m-n)]  can  be  exploited  in  inverting  the 
covariance  matrix. 

The  output  of  the  main  antenna  is  a  function  of  the  jamming 
signal  f(t)  and  of  the  scattering  process.  Representing  the  scattering 
by  a  large  set  of  independent  scatterers,  each  with  a  complex  reflection 
coefficient  a.  and  a  delay  t- , 

J  J 


y(t)  =  ^2  ctjf (t-Tj) 

j=i 

* 

The  corresponding  components  of  the  Ey  vector  of  (2)  are 


E„y*  -  En(t)y*(t) 


aj  p(Tj-nA) 

j-1 
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For  a  given  set  of  scatterers  represented  by  {o.,t.},  the  power 

J  J 

in  the  main  antenna  output  before  cancellation  is 


|y|2  ■  X)  aj  “k  p(VTj> 


This  result  follows  from  (4)  and  (6).  Note  that  this  is  also  the  first 
term  in  the  equation  (3)  for  the  output  power  after  cancellation. 
Combining  (3),  (7),  and  (8),  the  power  residue  in  the  output  is 


X)  vI!  [p(Tk‘Tj}  -  sj M_1  sk 

j,k=l 


i  M  S. 


where  S.  is  a  column  vector  associated  with  the  j  scatterer, 

J 


p(xj-A) 


sj  ■  p(Tr2A) 


p(Tj-NA) 


The  Monte  Carlo  technique  of  selecting  the  scattering  field  randomly 

and  running  the  program  many  times  would  yield  an  estimate  of  performance, 

using  (8)  and  (9)  for  the  power  before  and  after  cancellation,  respectively. 


A  more  efficient  method,  however,  is  to  assume  that  the  scatterers 


are  equally  spaced  with  independent  random  coefficients  {a-}. 

_  J 

ajak  =  c»  3  rea^  constant  for  j=k  and  0  for  j^k.  In  this  case 


|y|2  =  CJp(O) 


and 


E  [p(0) 

j=l 


*  -1 
si M  sj 


J 


The  cancellation  ratio  is 


C.R. 


=  1  - 


1 


J 


35157^  H'’sj 

j=l 


=  1 


N 


'mn  tnn 

m,n=l 


where 


The 


Hmn(j)  =  p(Tj-mA)p(Tj-nA)  j p(0 ) 
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3.  BANDL IMITED  NOISE  WITH  CONSTANT  SPECTRAL  DENSITY  OVER  THE  BAND^' 

The  performance  of  an  adaptive  scatter  canceller  depends  on  the 
spectrum  of  the  scattered  jamming  signal.  This  spectrum  depends  on  both 
the  spectrum  radiated  by  the  jammer  and  the  frequency  response  of  the 
receivers  in  the  receiving  system.  The  simplest  case  to  analyze  is  the 
bandlimited  rectangular  spectrum,  i.e.,  where  the  jammer  signal  is  strictly 
limited  to  the  bandwidth  B  with  a  constant  spectral  density  over  the  band. 
The  tap  spacing  is  assumed  to  be  A  =  1/B.  In  this  case,  the  autocorrelation 
function  of  the  jamming  signal  is 


P, (t)  =  f(t)f*(t-T)  *  0,(0)  05) 

The  components  of  the  E  vector  are  delayed  replicas  of  this  jammer 
signal,  independent  from  tap  to  tap  of  the  delay  line,  so  the  covariance 
matrix  M  is  diagonal  and  can  be  inverted  analytically.  In  this  strictly 
bandlimited  case,  the  residue  after  cancellation  is 
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It  was  shown  in  an  earlier  report  [1]  that  the  residue  of  (16)  goes  to 
zero  as  the  number  of  taps  N  becomes  very  large  in  the  bandlimited 
case  where  ( t )  is  given  by  (15).  However,  it  has  also  been  shown  [2] 
that  a  large  number  of  taps  extending  beyond  the  scatter  delay  interval 
is  required  in  this  bandlimited  case  to  achieve  a  large  cancellation 
ratio.  This  is  illustrated  in  Table  1.  For  example,  when  the  scattering 
region  extends  over  two  resolution  cells,  a  total  of  79  taps  yields  only 
25.9  dB  of  cancellation. 

Table  1.  Variation  of  Cancellation  Ratio  with  Number 
of  Taps  -  Bandlimited  Noise  with  Rectangular 
Frequency  Spectrum 


Number 

of 

Width  of  Scattering  Region  (2L) 

Taps 

(2N+1 ) 

2 

6 

14 

26 

38 

50 

70 

11.3 

2.9 

1.0 

0.5 

0.4 

0.3 

0.2 

15.3 

13.9 

2.9 

1.3 

0.9 

■29 

18.7 

18.5 

16.4 

3.7 

2.2 

1.5 

HIV 

21. 2 

21.2 

20.8 

18.3 

5.3 

3.3 

2.1 

23.3 

23.2 

23.1 

22.6 

21.3 

8.4 

4.1 

59 

24.6 

24.6 

24.6 

24.3 

23.9 

23.0 

7.9 

79 

25.9 

25.9 

25.9  |  25.7 

25.5 

25.2 

23.9 

4. 


GAUSSIAN  SPECTRUM  L 


Another  spectrum  which  is  often  assumed  in  analyzing  systems  is  the 
Gaussian  spectrum,  where  the  noise  spectral  density  has  the  form 

S(f)  =  e'6f2  (i; 

The  spectral  density  has  a  value  of  1/2  of  its  peak  value  of  unity  at 
f  =  .8325//B  .  The  3  dB  bandwidth  is 

B  =  (1j 

n 

The  auto  correlation  function  of  the  noise  process  corresponding 
to  the  Gaussian  spectrum  (17)  is 


p(t)  =  2  \  S(f )cos2Trfxdf 


■  /f 


1  p"Ut)VB 


(19 


If  the  delay  line  taps  (Fig.  1)  are  spaced  by  the  reciprocal  bandwidth, 

A  =  1/B  =  /b/1.665.  The  performance  of  the  scatter  canceller  has  been 

calculated  fora  variety  of  different  tap  spacings.  The  parameter  a  is 
the  ratio  of  assumed  tap  spacing  to  the  reciprocal  bandwidth,  i.e.. 


For  a  given  a  and  tap  spacing,  the  elements  of  the  covariance  matrix 

M  are 

Mmn  =  ^  exp  {-  [1T(m-n)Qt/1.665]2}  ,  (21 

where  the  parameter  8  has  been  assumed  equal  to  unity.  This  expression 

was  used  to  obtain  the  covariance  matrix  of  (5). 

The  elements  of  the  column  vector  S.  of  (10)  are 

J 

p(ij-nA)  =  /n  exp  | -[ir(Tj-na)/1.665)]^  |  ,  (22' 

This  equation  was  used  to  obtain  the  Hmn(j)  in  (14). 

Results  for  the  Gaussian  distribution  were  discussed  in  [4],  which 
contains  a  series  of  tables  showing  the  dependence  of  cancellation  ratio 
on  the  number  of  taps  and  the  tap  spacing.  One  set  of  results  is  shown 
in  Table  2,  illustrating  the  effect  of  the  number  of  delay  line  taps  on 
cancellation.  In  each  case,  the  tap  spacing  is  equal  to  the  width  of  a 
resolution  cell.  These  results  show  that  a  few  extra  taps  on  each  side 
of  the  scattering  interval  provides  good  cancellation.  With  the  scattering 
medium  extending  over  three  resolution  cells  of  delay,  a  total  of  7  to  9 
taps  yields  30  dB  of  cancellation. 


Table  2.  Effect  of  Number  of  Taps  on  Cancellation 

Gaussian  Spectrum 


N(Taps) 

L (Cells) 

K(Scat/Cell ) 

ot 

CR(dB) 

3 

3 

5 

.3333 

19.4 

5 

3 

5 

.3333 

27.1 

7 

3 

5 

.3333 

29.9 

9 

3 

5 

.3333 

31.0 

11 

3 

5 

.3333 

31.5 

Table  3  illustrates  the  effect  of  tap  spacing  on  cancellation. 
In  the  first  set  of  4  examples,  the  width  of  the  scattering  region  and 
the  interval  covered  by  the  taps  are  the  same,  viz,  an  interval  of  one 
reciprocal  bandwidth.  For  example,  with  a=.5  the  tap  spacing  and 
width  of  one  scatter  cell  are  both  1/2  of  the  reciprocal  bandwidth. 
Hence,  two  taps  and  two  cells  are  required  to  cover  one  delay  interval 
equal  to  a  reciprocal  bandwidth.  Note  in  this  series  of  4  examples 
that  the  cancellation  ratio  improves  significantly  as  a  decreases. 

The  next  set  of  4  cases  are  the  same  as  the  first  four  cases,  except 
that  two  additional  delay  line  taps  are  included. 
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These  results  show  that  small  tap  spacings  must  be  used  with  a 
Gaussian  spectrum,  e.g.,  a  £  1/3  to  achieve  ^  30  dB  of  cancellation.  In 
a  system  cancelling  scattered  interference  from  a  large  scattering  region, 
this  requirement  for  closely  spaced  taps  is  a  serious  problem.  Many 
extra  adaptive  degrees  of  freedom  would  be  required. 


Table  3.  Effect  of  Tap  Spacing  (ct)  on  Cancellation 
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5.  BANDL IMITED  COSINE  SPECTRUM 

Next,  consider  the  frequency  spectrum 


S(f )  =  ^  cos(^) 
=  0 


|f | <  B/2 
|  f  |  >  B/2 


(23) 


The  corresponding  autocorrelation  function  is 


R(t)  =  2 /S(f)  cos(2TtfT)  df 
o 

B/2 

=  J-  f  cos(^)  cos(2irfT)  df 


(24) 


From  (24),  it  can  easily  be  shown  that 


R(x)  =  cosJjtBt)  (25) 

1-(2Bt)^ 

Since  this  spectrum  is  strictly  limited  to  a  bandwidth  B,  it  suffices 
to  take  samples  at  an  interval  A  =  1/B.  With  samples  spaced  1/B  apart,  the 
exact  waveform  at  points  between  the  sample  points  can  be  reconstructed 
using  the  sampling  theorem.  The  elements  of  the  corresponding  covariance 
matrix  (5)  are 


The  column  vector  S,  for  the  jth  scatterer  is 


R(tj-A) 

R(xj-2A) 

• 

R(x.-NA) 

J 


(27) 


Again,  it  is  assumed  that  the  delay  line  taps  bracket  the  delay  interval 
of  the  scattered  jamming. 

The  performance  of  the  scatter  canceller  with  a  cosine  spectrum  was 
investigated  using  the  FORTRAN  program  listed  in  Appendix  1  of  this  report. 
The  cancellation  ratio  of  (13)  was  computed  for  a  set  of  different  cases. 

In  (13), the  covariance  matrix  M  was  obtained  from  (26),  and  Hmn  from 


Hmn(j)  =  R(xj-mA)  R(xj-nA)/R(o) ,  (28) 

where  R(x)  is  given  by  (25). 

First,  a  series  of  cases  were  run  to  determine  the  number  of  scatterers 


per  resolution  cell  required  to  simulate  a  continuously  distributed  scatter¬ 
ing  medium.  It  was  found  that  3  to  5  scatterers  per  cell  are  sufficient. 

A  series  of  examples  are  shown  in  Table  3  to  illustrate  this  point. 


Table  3.  Effect  of  Number  of  Scatterers  per  Cell  on  Cancellation  Ratio 


No.  Taps 

No.  Cells 

Scatterers/Cell 

CR/dB) 

6 

5 

1 

17.1 

6 

5 

2 

19.9 

6 

5 

3 

20.0 

6 

5 

4 

20.0 

6 

5 

5 

20.0 

6 

5 

10 

20.0 

8 

5 

1 

20.5 

8 

5 

2 

23.4 

8 

5 

3 

23.5 

8 

5 

4 

23.5 

8 

5 

5 

23.5 

8 

5 

10 

23.5 

The  correlation  function  for  a  cosine  distribution  (25)  goes 
to  zero  for  large  x  more  quickly  than  the  corresponding  sine  function  for 
the  rectangular  spectrum  (15).  One  would  expect  fewer  taps  bracketing  the 
scatterering  region  to  suffice  for  the  cosine  distribution.  The  Gaussian 
correlation  function  (19)  approaches  zero  more  rapidly  than  the  correlation 
function  of  (25),  so  more  bracketing  taps  should  be  required  with  the 
cosine  spectrum  than  with  the  Gaussian  spectrum.  These  expected  results 
were  confirmed  by  the  computer  analyses  of  the  3  cases. 

Table  4  shows  the  cancellation  ratio  as  a  function  of  the  number 
of  taps  for  the  case  where  all  the  scatterers  are  within  a  single  delay  cell 
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Table  4.  Effect  of  Number  of  Taps  on  CR  with  One  Scatterer  Cell 


No.  Taps 

No.  Cells 

Scatterers/Cel 1 

CR/  (dB ) 

2 

1 

5 

14.1 

4 

1 

5 

19.1 

6 

1 

5 

22.3 

8 

1 

5 

24.6 

10 

1 

5 

26.5 

12 

1 

5 

27.9 

14 

1 

5 

29.2 

16 

1 

5 

30.4 

18 

1 

5 

31.4 

20 

1 

5 

32.3 

and  the  cosine  spectrum.  Between  4  and  6  taps  are  required  for  20  dB  of 
cancellation,  and  14  to  16  taps  for  30  dB  cancellation.  Approximately  7 
taps  on  each  side  of  the  scattering  region  are  required  to  achieve  30  dB 
of  cancellation  in  the  case  of  Table  1. 

When  the  scattering  region  is  extended  over  several  resolution 
cells,  fewer  extra  cells  on  each  side  of  the  scatter  delay  interval  are 
required.  This  is  illustrated  in  Table  5  for  the  cases  of  4  and  5  cells 
of  scatterers.  In  these  cases  5  to  6  extra  taps  on  each  side  of  the  scatter 
interval  yield  30  dB  of  cancellation.  This  difference  from  the  results  of 
Table  1  is  expected, si  nee  with  more  scatter  delay  cells,  the  scattering  from 
the  central  cells  is  cancelled  more  exactly.  For  these  central  cells  there 


Table  5.  Effect  of  Number  of  Taps  on  CR 


No.  Taps 

No.  Cells 

Scatterers/Cell 

CR(dB) 

6 

5 

3 

20.0 

8 

5 

3 

23.5 

10 

5 

3 

25.7 

12 

5 

3 

27.5 

14 

5 

3 

28.9 

16 

5 

3 

30.1 

18 

5 

3 

31.2 

4 

4 

3 

15.1 

6 

4 

3 

20.9 

8 

4 

3 

23.9 

10 

4 

3 

26.0 

14 

4 

3 

26.0 

18 

4 

3 

31.2 

22 

4 

3 

33.0 

are,  in  effect,  more  taps  on  each  side  of  the  cells  which  can  be  used  to 
replicate  the  scattered  jamming  signal.  This  point  is  further  illustrated 
in  Table  6,  where  10  taps  yield  19.6  dB  of  cancellation  when  the  scatter 
delay  extends  over  10  cells.  With  8  cells  of  scatter  delay,  4  to  5  overlap 
cells  on  each  side  of  the  scatter  region  yields  30  dB  of  cancellation. 
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Table  6.  Results  for  Larger  Number  of  Scattering  Cells 


No.  Taps 

No.  Cells 

Scatterers/Cell 

CR(dB) 

10 

10 

3 

19.6 

12 

10 

3 

24.9 

14 

10 

3 

27.3 

9 

8 

3 

22.1 

11 

8 

3 

25.3 

13 

8 

3 

22.4 

15 

8 

3 

28.9 

17 

8 

3 

30.2 

6.  CONCLUSIONS 


In  an  adaptive  scatter  canceller,  the  achievable  performance  depends 
on  the  frequency  spectrum  of  the  signals  at  the  canceller,  the  number  of 
delay  taps  used,  and  the  spacing  between  taps.  The  complexity  of  the 
adaptive  system  depends  on  the  number  of  delay  taps,  which  may  be  large 
in  some  cases  of  interest  due  to  the  range  extent  of  the  scattering  medium. 
This  frequency  spectrum  is  determined  by  the  frequency  response  of  the 
receivers  in  the  usual  case  of  broad  band  noise  jamming,  and  the  receiver 
frequency  response  can  be  selected  to  facilitate  the  implementation  of  a 
scatter  canceller. 

Three  different  frequency  spectra  have  been  investigated,  i.e.,  the 
achievable  cancellation  calculated  as  a  function  of  tap  spacing  and  the 
number  of  taps  overlapping  the  scatter  delay  interval.  It  was  shown  that 
a  large  number  of  overlapping  taps  are  required  with  a  rectangular  band- 
limited  spectrum.  With  a  Gaussian  spectrum,  the  required  tap  spacing  is 
roughly  1/3  of  the  reciprocal  bandwidth,  thus  requiring  a  large  number  of 
taps.  With  a  cosine  spectrum,  a  tap  spacing  of  the  reciprocal  bandwidth 
can  be  used  and  a  few  overlapping  taps  yield  20  to  30  dB  of  scatter  can¬ 
cellation.  The  cosine  spectrum  will  be  assumed  in  continuing  studies  of 
the  adaptive  scatter  canceller.  These  studies  will  address  the  important 
problem  of  designing  a  practicable  adaptive  weight  computer  in  systems 
where  the  scatter  delay  interval  extends  over  a  large  number  of  resolution 
cells. 

In  the  first  quarterly  report,  it  was  shown  that  the  adaptive  scatter 
canceller  can  be  used  when  two  jammers  are  illuminating  the  scattering  region 
in  the  main  beam.  This  requires  the  use  of  two  auxiliary  antennas  which 
receive  the  direct  path  jamming  signals^. 
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APPENDIX  1  .  SCATTER  CANCELLATION  WITH  COSINE  SPECTRUM 


FORTRAN  IV  H01A-1  THU  05-AUG-82  00:00:00 


0001 

C 

C 

C 

0002. 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

c 

c 

C 

0010 

0011 

0012 

0013 

0015 


0018 

0019 


0020  110 
c 
c 
c 

0021 

c 

c 

C 

0022 

0023 

0024 

0025 

0026 

0027 

0023 

0029 

0030 

0031  130 

0032  D 
0033  120 


PROGRAM  MAIN 
NN-NUMBER  OF  TAPS 
LL-NUMBER  OF  SCATTER  CELLS 
KL-NUMBER  OF  SCATTERERS  PER  CELL 
DIMENSION  B 1 50, 50) 

DOUBLE  PRECISION  A I 30 , 30 ) , H I 30 , 30 ) , CR , T AU , A 1 , A2 
WRITE (7, 100) 

ACCEPT  200 , NN , LL , KL 
PI-4.*ATAN< 1. ) 

NM- NN 

DO  150  1 1  - 1  , NM 
NN=LL+2*I 1 1  — 1 ) 

FORM  COVARIANCE  MATRIX 


DO  110  M-1,NN 
AIM,M)*=1. 

BIM,M)-AIM,M) 

IF(M.EQ.NN)  GO  TO  110 
DO  110  N-M+l , NN 

AI M,  N) -COSIPI* I M-N) ) / ( 1 . -4. *IM-N) **2) 

AIN, M)=AIM,N) 

BIM, N)=AIM,N) 

BIN, M)=AIN,M) 

CONTINUE 

INVERT  COVARIANCE  MATRIX 

CALL  MATRIA,NN, 1,30) 

COMPUTE  H(M,N) 

Tl-FLOAT (NN-LL+1 ) /2. 

N1~LL*KL 
DO*  120  M- 1 ,  NN 
DO  120  N-l , NN 
HIM, N)-0. 

DO  130  J-l ,  N1 

TAU=T 1 + ( FLOAT I J ) - . 5 ) /KL 

A 1 -DCOS ( P I *  <  T AU-M ) ) / ( 1 . -4. *<TAU-M>**2> 

A2=DC0S ( P I * ( TAU-N )  ) / 1 1 . -4 . *  I  TAU-N ) **2 ) 

HIM, N) -HIM, N) +A 1 *A2 

WRITE (6, 800)  M, N, HIM, N) , AIM, N) , B(M, N) 
CONTINUE 


COMPUTE  CANCELLATION  RATIO 


c 

c 


0034 

C2- 1 . /FLOAT (N1 ) 

0O35 

CR-1 . 

0036 

DO  140  M= 1 . NN 

0037 

DO  140  N-1,NN 

•:*033 

140  CR=CR-C2»A(M,N) »H(M, N) 

OkO'? 

IF  (  CR.  LT.  0.  )  CR=.  0 . 00 -»1 

O  ■'  •’  1 

CRDT-  =- 1 0'-TLOG  1  <'» (  CR ) 

tvl 


i 

0042 

0043 

0044 

150 

iS 

0045 

100 

>• 

0046 

200 

0047 

500 

•  ■ 

004© 

700 

004? 

SCO 

OO50 

WRITE (6, 500  >  NN, LL, KL 
WRITE <6, 700 )CRDB 
CONTINUE 

FORMAT </X, 'INPUT  #  TAPS,#SCAT  CELLS*  #SCAT /CELL  IN  315' 
FORMAT (315) 

FORMAT < /X, '#TAPS=' , 14, 5X, ' &CELL3=' , 14, 5X, ' #SCAT/GELL= ‘ 
FORMAT (/X,  ••CANCELLATION  RATIO (DB)='  ,F10.4) 

FORMAT ( X , 2 1 1 0 » 1 ©X , 3F 1 5.  5 ) 

END 


|  * 


fc  K‘> 


mi 


,14) 


